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Abstract 

Color fluctuations in hadron-hadron collisions are responsible for the presence 
of inelastic diffraction and lead to distinctive differences between the Gribov 
picture of high energy scattering and the low energy Glauber picture. We find 
that color fluctuations give a larger contribution to the fluctuations of the num- 
ber of wounded nucleons than the fluctuations of the number of nucleons at a 
given impact parameter. The two contributions for the impact parameter aver- 
aged fluctuations are comparable. As a result, standard procedures for selecting 
peripheral (central) collisions lead to selection of configurations in the projectile 
which interact with smaller (larger) than average strength. We suggest that 
studies of pA collisions with a hard trigger may allow to observe effects of color 
fluctuations. 



1. Introduction 

Currently most of the experimental studies as well as modeling of the nucleus- 
nucleus (proton-nucleus) collisions involve using the Glauber model. Namely, 
the number of involved nucleons is calculated probabilistically assuming that 
each Nucleon-Nucleon (NN) inelastic collision is determined by the value of 
a^ n N at the collision energy. 

However, the dominance of large longitudinal distances in high energy scat- 
tering [l[ changes qualitatively the pattern of multiple interactions. Indeed, 
in the Glauber approximation high energy interactions of the projectile with a 
target occur via consecutive rescatterings of the projectile off the constituents 
of the target. The projectile during the interactions is on mass shell - one takes 
the residues in the propagators of the projectile. This approximation contra- 
dicts the QCD based space-time evolution of high energy processes dominated 
by particles production. The projectile interacts with the target in frozen con- 
figurations since the life time of the configurations becomes much larger than 
the size of the target. Hence there is no time for a frozen configuration in the 
projectile to combine back into the projectile during the time of the order Rt, 
the radius of the target. As a result the amplitudes described by Glauber model 
diagrams die out at large energies a 1/s (a formal proof which is based on the 
analytic properties of the Feynman diagrams was given in [H, HJ). 
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In the Glauber model the number of interacting nucleons is calculated prob- 
abilistically assuming that the probability of individual NN inelastic collisions 
is determined by the value of cr^ n N at the collision energy. Fluctuations of the 
number of wounded nucleons at a given impact parameter are generated solely 
by fluctuations of the positions of nucleons in the nucleus and (in some models) 
due to peripheral collisions of nucleons, where the interaction is gray and hence 
the chance to interact differs from one or zero. Hard collisions are treated as 
binary collisions, which is equivalent to taking the diagonal generalized parton 
densities of nuclei, Ja(x, Q 2 ,b), proportional to the impact factor T(b): 

F A (x,Q 2 ,b) = f N (x,Q 2 )T(b), (1) 

where T(b) is normalized as J dbT(b) — A. A nuclear shadowing correction is 
introduced for x < 0.01. 

The high energy theory of soft interactions with nuclei was developed by 
Gribov Q who expressed the shadowing contribution to the cross section of 
hadron-nucleus (hA) interactions through the contribution of non-planar dia- 
grams. The Gribov- Glauber theory, in difference from the low energy Glauber 
theory, requires taking into account that a particular quark-gluon configuration 
of the projectile is frozen during the collision and that it may interact with dif- 
ferent strength as compared to the average strength. This leads to fluctuations 
of the number of collisions which are significantly larger than in the Glauber 
model. The fluctuations of the strength of the interaction are related to the 
ratio of inelastic and elastic diffraction in NN scattering at t = 0. Relevance of 
fluctuations of the strength was first pointed out in Q but these effects were 
never analyzed in detail before. 

Another effect contributing to fluctuations of observables in hA collisions is 
fluctuations of the gluon density which can originate both from the fluctuations 
of the nucleon configurations and from the fluctuations of the gluon densities in 
the individual nucleons. We will consider this effect elsewhere. 

The paper is organized as follows. In section [2] we summarize the necessary 
information about fluctuations of the strength of NN interaction. In section[3]wc 
use the Gribov- Glauber model in the optical approximation to obtain analytic 
results for the strength of fluctuations of the number of wounded nucleons and 
relative contributions to these fluctuations of the fluctuations of the strength 
of the interaction and of geometry of collisions. In section [4] we develop a full 
Monte Carlo (MC) model in which the geometry of projectile-target nucleon 
interaction is accounted for, and the strength of the interaction fluctuates on 
an event-by-event basis. The results for the change of the distribution over 
the number of collisions and for the dependence of the average strength of the 
interaction on impact parameter are presented. Possibilities for observing color 
fluctuation effects in collisions with hard triggers are outlined. 
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2. Color fluctuation effects in proton-nucleus collisions 

2.1. Gribov inelastic shadowing 

It was demonstrated by Gribov [3] that the nuclear shadowing contribution 
to the total cross section of the hadron-deuteron scattering can be expressed 
through the diffraction cross section at t = 0. Operationally this amounts 
to the replacement in the Glauber formulae of the elastic hN cross section at 
t ~ by the sum of elastic and diffractive cross section at t = 0, leading to 
an enhancement of the multinucleon interactions. For heavier nuclei the Gribov 
formulae involve the coupling of the projectile to N > 2 vacuum exchanges which 
has to be modeled. 

The contribution of the double scattering to the total hadron-nucleon (hN) 
cross section is enhanced by a factor 1 + u> a , where 



The relation between the double scattering cross section and the total diffraction 
cross section can be naturally understood in the Good and Walker formalism [A, 
which provided the effective realization of the Fcinbcrg-Pomeranchuk picture [8| 
of the inelastic diffraction. In this formalism one introduces cigenstates of the 
scattering matrix diagonal in a; see Ref. @ for a review. Configurations with 
different o~i scatter without interference off two target nucleons contributing in 
the case of scattering of two nucleons with strength oc of to the shadowing of 
the total cross section. This is the same quantity as in the expression for the 
total cross section of hadron-nucleon diffraction at t = 0. This interpretation 
of the Gribov result for the shadowing correction to the total cross section was 
first given by Kopeliovich and Lapidus (loj . 

2.2. Distribution over the strength of interaction 

The fluctuations of strength of interaction arise naturally in QCD where the 
strength of interaction depends on the volume occupied by color. In particular, 
the presence of some small configurations leads to fluctuations interacting with 
a small cross section. So we will refer to these fluctuations as color fluctuations. 

In order to describe the effect of color fluctuations for a variety of processes 
it is convenient to introduce the notion of distribution over the strength of 
interaction, Ph(o~tot) - the probability for an incoming hadron to interact with 
total cross section 17404. The distribution Ph{o~tot) satisfies two normalization 
sum rules: 







and the Miettenen-Pumplin relation 11 




(4) 
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where is the free cross section. Experimentally, w CT first grows with energy 
then starts dropping at energies ^/s > 100 GeV. There are no direct measure- 
ments at the RHIC energy of 200 GeV, but an overall analysis indicates that it 
is of the order 0.25. The first LHC data seem to indicate that inelastic diffrac- 
tion still constitutes a large fraction of the cross section - it is comparable to 
the elastic cross section, suggesting uj a ~ 0.2 at those energies. It is difficult 
at the moment to ascribe error bars to these numbers. However, it is expected 
that the values of u a corresponding to the LHC energies will be soon measured 
with a good precision. 

It is worth emphasizing here that these seemingly small values of cu a corre- 
spond to very large fluctuations of the interaction strength. For example, if we 
consider a simple two component model (equivalent to the quasi-eikonal approx- 
imation), in which two components are present in the projectile wave function 
with equal probability and interact with strengths <j^ t and a^ t : 

ail] = - v/SJ) , ag> = o*£(l + v/£J) . (5) 

Thus for u a = 0.25, we have crtot/ a tot = 0.5, er^/cr^ = 1.5 and hence 

" tot I "tot — °- 



3. Gribov - Glauber model predictions for fluctuations in the optical 
approximation 

In order to illustrate the effects of the color fluctuations and their interplay 
with the fluctuations of the local nuclear density we first consider the optical 
approximation of the Glauber model where the radius of the NN interaction is 
neglected as compared to the distance between the nucleons. 

Within this model the total inelastic hadron-nucleus cross section can 
be written as follows: 

= / * (i - [i - m A ) - E J *>*<® N . ( 6 ) 



where x(b) = of T f T(b)/A and normalization / dbT(b) = A. 

Note that in Eq. (|6]) nucleon-nucleon correlations in the nuclear wave func- 
tion arc neglected as well as the finite radius of the hadron-nucleon interaction; 
an implementation of correlations in the optical limit in the Gribov - Glauber 



formalism can be found in Rcfs. [12|, [13j and in Ref. [14[ within the MC ap- 



proach and will not be discussed here. Eq. © can be rewritten as a sum of 



positive cross sections [15j as follows 



°l A = E ^' ^ = {A - A N)IN\ J dbx(b)N [1 " x{b)]A ~ N ' (7) 

where gm denotes the cross section of the physical process in which N nucleons 
have been involved in inelastic interactions with the projectile. Using Eq. (J7J), 
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the average number of interactions (N) can be expressed as 



A I A n hN P Arr hN 

(N) = £ AW = / dbT(b) = * , (8) 

jV=l ' N=l in •> a in 

which coincides with the naive estimate of shadowing as being equal to the 
number of nucleons shadowed in a typical hA inelastic collision. 

Wc can include color fluctuations by allowing the inelastic cross section a in 
to be distributed according to a proper distribution, P H (<r in ) : 

u hA = J da m P H (a m ) J db (l - [1 - x{b)] A ^ , (9) 

where now x(b) — Oi n T(b)/A, and 

°N = j d ^ P H(^n) {A jl y Nl J dbx(b) N [1 - X (b)] A ~ N . (10) 

The probability of collisions with exactly N inelastic interactions in both Glauber 
model and the color fluctuation approximation is simply Rn = /&i A - 

Using the equations above we can for example calculate the average number 
of collisions which is given by the same equation as for the Glauber model (Eq. 
@), leading to a very small (a few %) change of average N, as shown in Table 
[U since the inelastic corrections to tr^ 4 are small for a realistic Piii&in); see 
Ref. [lH and references therein. The physical reason why the corrections are 
small is that, in a broad range of b, the interaction is close to the black limit 
for all essential values of ai n , so only a small range of (large) b contributes to 
inelastic shadowing corrections. At the same time the color fluctuation effect is 
large for the variance of the distribution over the number of collisions. Eq. (|10[) 
leads to 



(N(N-l))=A(A-l)^ I dbT\b), 
and hence the variance is equal to 

^^-! = ^^i# ( d bT\b) + '-l. (12) 
(N) (N) <rt A J W 



One can see from Eq. (jT2]) that the variance receives contributions both from the 
fluctuations of the impact parameter and from the fluctuations of ai n . Using 
Eqs.(|8]),([II]) we obtain for the variance in Eq.(fT2]) the value of about 0.46 (RHIC) 
and 0.51 (LHC). Numerical values of the different terms in Eq.lfTS"]) are: 1.26 
+0.20 -1 = 0.46 (RHIC) and 1.38 +0.13 -1 = 0.51 (LHC). The account of 
the color fluctuations practically does not change (N). It mainly changes the 
nominator of the first term by the factor l + w CT o Though this change is rather 



x We assume here that fluctuations for the inelastic and total cross sections are similar, cf. 
discussion before Eq. 11171 . 



5 



small, the strong cancellation between the first and the third terms of Eq. (|12p 
strongly enhances the effect of color fluctuations. 

A more realistic treatment of the color fluctuations taking into account the 
profile function of the NN interactions and small effect of short-range corre- 
lations is possible in the MC model described in the next section. First, one 
calculates the probability Pjv(&) shown in Fig. [T]of having exactly A inelastic in- 
teractions at a given impact parameter b. Next one can calculate the quantity in 
Eq. (fT2"| by integrating Pm (b) over the impact parameter: Pjy = 27r J b db Pm(b). 
The results are given in Table [TJ 





Monte Carlo 


Optical Model 


energy /model 


(AO 




U1 N 


(A) 






RHIC, Glauber 


4.6 


31.6 


0.51 


5.0 


35.9 


0.46 


RHIC, GG2 


4.7 


38.9 


0.74 


5.1 


45.3 


0.71 


RHIC, GG P h {atot) 


4.8 


39.2 


0.72 


5.2 


45.6 


0.70 


LHC, Glauber 


6.7 


72.4 


0.59 


7.6 


88.0 


0.51 


LHC, GG2 


6.8 


84.2 


0.80 


7.8 


106.2 


0.75 


LHC, GG P h {a tot ) 


6.8 


82.1 


0.77 


7.8 


106.4 


0.74 



Table 1: The fluctuations, as defined in Eq. 1121 . calculated both within the MC approach 
and optical model. We used no color fluctuation (Glauber), color fluctuations implemented 
with the two states model described in the text (GG2) and with the full color fluctuation 
model (GG Ph{&tot)) described by the distribution Ph(ctot) of Eq. J16t . 



A comparison of some of the predictions of the optical approximation of the 
Glauber model and the MC calculations, which take into account finite radius 
of the NN interaction neglected in the optical model, will be given below. 

4. Monte Carlo algorithm for modeling effects of fluctuations 

Wc have seen from the analysis of the optical model that fluctuations in the 
number of wounded nuclcons originate both from color fluctuations and from 
fluctuations of the number of nucleons along the path of the projectile. 

The event-by-event fluctuations of the number of wounded nucleons due to 
the fluctuations in the number of nucleons at a given impact parameter are 



present already on the level of the Glauber model [141 ]. In the case when no 
fluctuations of a are present, (A(cr l ' l ^ v )) is given by Eq. ((HJ). In this case we can 
write 

(N(a™) 2 } = {N) 2 (l+u p (a™)), (13) 

where co p (a!^) is the dispersion in the case of no color fluctuations. We found 
that uip(a^) drops as a function of cr^ , as a consequence of the increasing num- 
ber of nucleons in the interaction volume. In the calculations we use the event 



generator [14j . This event generator includes short-range correlations between 
nucleons, however this effect leads to a very small correction for the discussed 
quantity. The code also includes a realistic dependence of the probability of 
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the NN interaction on the relative impact parameter of the projectile b, and 
the target nucleon b 3 ■: b — bj . The probability of the interaction is expressed 
through the impact factor of the NN elastic amplitude 

r(6 - 6,) = e -(b-b 3 f/2B (M) 

as follows: 

p(b,b,) = i - [i - r(6-b,-)] 2 • (is) 

Here we used the exponential fit to the elastic cross section da/dt oc exp(Bt). 

In order to perform numerical analyses we follow (l6| . and take the proba- 
bility distribution for cr tot as follows: 

D i \ / °~totl°~<d - 1\ , 1g n 

Ph(<7tot) = P : exp<^ -2 > , (16) 

where p is a normalization constant and we have Co =72.5 mb and Q =1.01 
at LHC energies, while <tq =32.6 mb and fl =1.49 at RHIC energies. One can 
verify that the distribution of Eq. (|T5)) satisfies the sum rules (j3J), ©, with our 
values o*£ = ag? = 51.95 mb for RHIC and a[f t = o*™ = 94.8 mb for LHC 
energies. 

When converting from the distribution over a tot, Ph{vtot), to the distribution 
over am, Pni^in), we used the geometric scaling observation that the t-slope of 
the elastic scattering is proportional to a tot- So the ratio ain/ctot — A weakly 
depends on the projectile and energy. Hence we take A = const, so that we 
simply have to use a Jacobian 1/A, with 

PH((?in) = Ph(crtot)/^ , cr in = A a to t ■ (17) 

Indeed in this case J do~i n Pn {<Ti n )~\ holds as well. This corresponds to B{a to t) = 
B K a tot ) °tot/o- t ot ■ 

In our numerical studies we used the fluctuation distribution given by Eq. 
C5K, Otot* S iven above and B=U GeV~ 2 (RHIC), B = 19.38 GeV" 2 (LHC). 
This parametrization satisfies the s-channcl unitarity condition L(6) < 1. In our 
model this condition holds automatically also for the elastic "color- fluctuation" - 
nucleon amplitude. Our algorithm is a natural extension of that of [l4^ - where 
distribution over N was calculated in the Glauber model neglecting effects of 
color fluctuations. 

Since the contributions of states with different a do not interfere, the prob- 
ability Pn(P) to have exactly N inelastic interactions at given b ifH 

P N (b) = / dat tP h (a tot )P N (b;atot), (18) 



2 In this treatment we neglect small contributions of incoherent diffractive processes pA — * 
XA* , which mostly contribute to Pi(b). 
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where P/v(&;<7tot) is calculated using the procedure of Ref. [14[ for fixed 
in the Glauber model. Including color fluctuations results in a substantially 
broader distribution over b of the probability P/v (b) of having exactly N inter- 
actions for a given impact parameter N, as shown in Fig. [1] The two-component 
model gives the distributions pretty close to the distributions including full fluc- 
tuations. P/v (6) are obviously normalized so that J2 N J G?6P/v(fr) = a^. The 
calculations of Table [T] have been performed integrating the quantities of Fig. 
[1] over the impact parameter: P/v = / dbP/v(6); (TV) = J2n ^Pn/ ^2 N Pn', 

(n 2 ) = E n n 2 Pn/E n Pn. 

Another quantity which characterizes the effects of spatial and color fluctu- 
ations is dispersion of the number of interactions at a given impact parameter, 
b. To illustrate the expected pattern let us first consider the case of small b 
and large A, when the probability of having at least one inelastic interaction is 
1. In this case (N) = T(b)(Ti n , hence the dispersion of the distribution over N 
including both effects can be calculated as follows: 

(N 2 ) = J da in P H (a in ) (N) 2 f T^y) (1 + (<*»)) ■ (19) 

Now we can calculate the total dispersion. The first term in (1 + u) p ) simply 
gives Lu a . The second term takes into account the dependence of oj p on the 
fluctuating <Ji n : 



J da in P H (crin) fz^ry) uj p (<r in ) . (20) 



Since the integral in the second term is dominated by <Ji n > , for which ui p 
is smaller than in correspondence of the average value of <7j n , , Eq. (|20[) 
leads to a dispersion somewhat smaller that io a + cj p (a^). This is consistent 
with the pattern we find in the numerical calculation presented in Fig. [5] for 



(N 2 ) b - (N)l 



Dib) - vjv ib ,~r ib , (2D 



(N) b = J2 N NP N (b)/J2 N P N (b) and (N^) b = J2 N N 2 P N (b)/ J2 N J P N (b). One 
can see that for RHIC and LHC energies the dominant effect comes from color 
fluctuations. Moreover, the two states approximation gives the result which is 
very close to the calculation with full P/ l (ut t), so the two states model can be 
used to simplify modeling of color fluctuation effects. 

The large variance of the distribution leads to a much wider distribution 
over TV than in the Glauber model, as shown in Fig. |3] The figure shows the 
quantities P/v = / db P/v (b) / a^f ; the same quantities are plotted in logarithmic 
scale in the insets, and one can see that the color fluctuations produce a much 
stronger large N tail. Among other things, this implies that selection of events 
which in the Glauber model correspond to very central impact parameters ac- 
tually gets a significant contribution from pretty large impact parameters - for 
example, in the two component model discussed above the collisions at impact 
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LHC - N=3 




LHC - N=7 




LHC - N=10 

Glauber - 
Gribov-Glauber two states 
Gribov-Glauber P(a) - 




2 4 6 

b[fm] 

Figure 1 : The probability Pjv (&) of having N inelastically interacting (wounded) nuclcons in 
a pA collision, vs. impact parameter b, when using simple Glauber (red curves), a two states 
model (black curves) and a distribution P h (crtot) (blue curves); cf. Eq. H16I I. The Pj^(b)'s 
are obtained by extension of the MC code of Ref. Il4l to include color fluctuations. Top row 
shows Pjv=i(fe); the remaining panels correspond to N = (N) and N = (N) ±0.5{7V). (N) is 
taken as 5 and 7 for RHIC and LHC energies, respectively (cf. Table [!}■ 
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Figure 2: Effect on fluctuations of the dispersion, Eq. (!_' I li . when using a distribution of atot 
with two values of the cross section with equal probability and with Ph{o~tot) given by Eg. J161 I. 
for realistic parameters corresponding to RHIC (left) and LHC (right) energies. 



parameter b satisfying the condition T(b)/T(0) = 1/(1 + ^/uj) with a probability 
of 1/2 generates the same number of wounded nucleons as average number of 
wounded nucleons at b = 0. For us = 0.25 we have 1/(1 + s/uj) = 0.67 and this 
corresponds to b ~ 4.58 fm. 

An important implication of the broad distributions over N which is mostly 
due to fluctuations of the strength of the interaction is that selection of large N 
also selects configurations in the projectile nucleon with cross section larger than 
average. To illustrate this trend within our MC, let us consider the average a tot 
for events with a given number N of wounded nucleons. Denoting the probability 
to have exactly N wounded nucleons Pn = J dbPjsi(b) and using Eq. (fl"8|) . we 
can write 



0.25 
0.2 

0.15 
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0.05 



RHIC energy 
0.1 




5 15 25 

Glauber 

Gr.-Glauber two states 
Gribov-Glauber P(e>) — 
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N (Wounded Nucleons) 
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Glauber 

Gr.-Glauber two states 
., Gribov-Glauber P(a) — 





10 15 20 

N (Wounded Nucleons) 



Figure 3: Effect of the event- by-event fluctuating values of atot, for RHIC (left panel) 
and LHC energies (right panel) on the number of wounded nucleons, calculated as = 
J db Pjv (b) I a 1 ?^ . Red curves show the results obtained with the usual Glauber calculation 
with fixed cross section, black curves correspond to calculations with the two-component atot 
model and blue curves correspond to calculations with fluctuating cross section with Ph{<?tot) 
distribution. The insets show the same quantities in logarithmic scale. 



(ctotjN 



1 j datotdb a to t Ph(a- to t)PN(b; atot) 



&tot I dvtotdb P h {a t at) Pw{b\ a t0 



(22) 
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The results of the calculation are presented in Fig. H) One can see that selecting 
N 3> (N) leads to a significant enhancement of the contribution of configura- 
tions which have interaction strength larger than average. For small N average 
(ctot) n is below aprf , but the effect is relatively small especially for N — 1 where 
very peripheral collisions contribute which are not sensitive to the fluctuations. 
A natural source of large a's are configurations of larger than average transverse 
size. One can expect that the gluon field is enhanced in these configurations 
while the distribution in x - the light-cone fraction carried by partons of the 
projectile - is softer for large x leading to a correlation between the distribution 
over N and distribution over a; of a hard collision. 



RHIC energy LHC energy 




0.9 



2 4 6 8 10 12 14 16 18 20 5 10 15 20 25 30 

N (Wounded Nuoleons) N (Wounded Nucleons) 



Figure 4: Effect of fluctuations on the event-by-event fluctuating values of atot, for RHIC and 
LHC energies. 



Matching the number of wounded nucleons to the physical obscrvables is 
certainly a challenging problem in view of fluctuations of the impact parameter 
in the collisions. A model independent treatment of this problem would require 
a study of pA collisions for different nuclei. Still the central multiplicity appears 
to be a good observable even in the presence of the color fluctuations. Indeed 
in the soft interaction dynamics the hadron multiplicity for central rapidities, 
Vc.m. ~ 0, does not depend on , as it is determined by the density of 
partons in a single Pomeron ladder. Hence the hadron multiplicity for y c . m . ~ 
should be about the same for different fluctuations. Also the first studies 
of the pA collisions at the LHC indicate that to a good approximation the 
hadron multiplicity for pt > 1 GeV is proportional to the number of wounded 



nucleons calculated in the Glauber model [17[. Hence we expect that selecting 
events with the y c . m . ~ hadron multiplicities: Mj (M) > 2.5 should select 
configurations in the projectile significantly larger than average ones (cf. Fig. 
Hp) with significantly different parton distributions. 

Correspondingly, a trigger for configurations of smaller than average size 
would lead to a more narrow distribution in TV. One such possibility is to select 
as a trigger a hard process in which a parton of the proton with x p > 0.6 is in- 
volved. One may expect that in this case one selects quark-gluon configurations 
without qq pairs and significantly screened gluon field, leading to a in signifi- 



cantly smaller than average and hence a strong suppression of large N tail [18 1 
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Such measurements appear to be feasible using the data collected in the 2013 
pA run at the LHC in which a significant number of events with large x p should 
have been collected. Since this kinematics (for the current LHC detectors) cor- 
responds to very large pr's of the jets, one expects that for the inclusive cross 
section impulse approximation would work very well. Hence it would be possible 
to avoid issues of the final / initial state interactions and nuclear shadowing in 
interpreting these data. 

A convenient quantity to study these effects experimentally would be a mea- 
surement of the distribution over x p for different classes of hard collisions at 
fixed xa normalized to the distribution in the inclusive pA scattering. A large 
effect is expected for the central collisions where the hard cross section should 
be suppressed for large x p > 0.2 ~ 0.3 and enhanced for x < 0.05. 

Note that such a measurement among other things would allow to test in an 
unambiguous way the explanation of the EMC effect at large x as due to the 
dominance of the smaller than average size configurations in nuclcon at x > 0.6; 
for a recent review see Ref. [lij . 

We also investigated the impact of fluctuations of the definition of ccntrality 
classes. We followed the experimental definition, in which the centrality is 
proportional to the fraction of total inelastic cross section provided by a given 
type of events. We can extract from the MC results of Fig. Q]thc probability Qm 
of having at least N inelastic interactions, irrespective of the impact parameter 
b (cf. Eq. 0): 

Qn ~ 7— , (23) 

J2m=i J dbP M {b) 

in such a way that Qn=i — 1 by definition. This allows to estimate the fraction 
of er^ 4 arising from a given interval in the number of wounded nucleons. Then, 
one can choose a centrality class and select the interval in number of wounded 
nucleons which contributes to that class. In Fig. [51 we have chosen the classes of 
the 20% most central events by requiring it to provide 20% of the total inelastic 
cross section and, similarly, we have singled out the 20%-40% and 40%-60% 
centrality classes, and the 40% most peripheral events as the last class. We use 
the number of the wounded nucleons corresponding to (closer to) these cuts as 
limits in N entering in Eq. ([24|) . for the calculation of the curves in Fig. El In 
Fig. [5] we show, for the selected classes, the distribution of events as a function 
of impact parameter by plotting 

bF(b) = b P ^( fo )> ( 24 ) 

N=N min 

where N m i n and N max are the values singled out by the cuts described above 
and shown in Fig. [5j The distribution of Fig. |6] was calculated both with the 
usual Glauber approach, i.e. with a fixed and with the inclusion of a 

fluctuating cross section according to Phip'tot)- For all the classes and both the 
RHIC and LHC energies, it can be seen that the fluctuating cross section tends 
to push the distribution toward larger values of the impact parameter. 
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Figure 5: Fraction of inelastic cross section plotted as a distribution over impact parameter 
as defined in Eq. (123 D - Horizontal lines at 0.2, 0.4 and 0.6 correspond to the experimental 
definition of 20%, 40% and 60% centrality, respectively. 



5. Conclusions 

Wc have demonstrated that color fluctuations lead to a significant modifica- 
tion of the distribution over the number of nucleons involved in inelastic proton 
- nucleus collisions at collider energies. Study of the correlations between the 
soft central multiplicity and the rate of hard parton-parton interactions in the 
pA collisions at the LHC would provide a new avenue for investigating the three 
dimensional structure of proton. In particular such measurement will allow to 
test a conjecture that quark-gluon configurations in the proton containing large 
x p partons have a significantly smaller than average size. 
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